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Abstract 

The rates of axion emission by nucleon-nucleon bremsstrahlung are calculated 
with the inclusion of the full momentum contribution from a nuclear one 
pion exchange (OPE) potential. The contributions of the neutron-neutron 
(nn), proton-proton (pp) and neutron-proton (np) processes in both the non- 
degenerate and degenerate limits are explicitly given. We find that the finite- 
momentum corrections to the emissivities are quantitatively significant for 
the non-degenerate regime and temperature-dependent, and should affect the 
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existing axion mass bounds. The trend of these nuclear effects is to diminish 
the emissivities. 

Key words: nuclear effects, bremsstrahlung, axions, neutron stars 
PACS: 14.80.Mz, 97.60. Jd, 95.30Cq 



1. Introduction 

The search for new particles/interactions beyond the Standard Model is 
one of the most important tasks of particle physics. While several candi- 
dates and proposals may be considered as "exotic" , in the sense of not being 
required by the data, it is generally agreed that there are some possibilities 
definitely expected as minimal extensions of it. Axions (El- jij]) belong to the 
latter category, as expected from the Peccei-Quinn ([3[]-[5|) solution to the 
strong CP-problem. Axions are pseudo Nambu-Goldstone boson associated 
with the spontaneous breaking of the Peccei-Quinn symmetry. Their masses 
and couplings are directly related to this symmetry-breaking scale. Viable 
versions of axionic models include the KSVZ 

(Q,H) and DFSZ (0-0) 

ax- 
ions coupled to hadrons only and to leptons and hadrons respectively. 

In addition to experimental efforts for a direct detection of axions, astro- 
physical and cosmological arguments have played a key role in their search. 
Actually, stringent bounds have been obtained from the consideration of 
horizontal branch stars (0- 11]), white dwarf cooling and SN1987A neutrino 



pulse duration f [12| N ). among others. A general review of these arguments 



has been given in [13| (see also [14j-[16j for a thorough account 



One of the main ingredients for an accurate calculation of axion mass 
bounds is the emissivity in the nucleon bremsstrahlung reaction NN — > 
NNa, thought to be dominant in important astrophysical events, such as 
newly born neutron stars. Calculations and discussions on the applicability 



of the emissivity formulae were given in [8j,[12j-[19j. Quite generally, in these 
papers, the calculations were performed for a one-pion exchange free nucleon 
gas, leaving aside important issues later incorporated and assessed, such as 
the effects of correlations between nucleons (see for instance [20]). Attempts 
to link the emissivity to laboratory data have been also made. We would like 
to present in this communication a reassessment of the bremsstrahlung emis- 
sivity including the full momentum dependence of the matrix elements. We 
found that the hitherto neglected dependence produces large temperature- 
dependent corrections to the rates independently of the many body effects, a 



2 



feature that points by itself to a revision of some of the astrophysical bounds. 

Axion emission are important for the evolution of stars particularly for 
(hot) neutron stars (NS). For the conditions relevant to the core of hot NS just 
after their formation (T ~ 30 — 60 MeV, p > po (with po = 2.7 x 10 u gcm~ 3 
the nuclear matter density) the dominant emission process are the nucleon- 
nucleon (NN) bremsstrahlung ([191]) 



n + n — > n + n + a (a) 
p + p — > p + p + a (b) 
n + p^n + p + a (c) 
Previous calculations of the axion emissivities have been performed by 



Iwamoto ([9j,[l7j) in the degenerate (D) limit. For the NN interaction he 
used a OPE potential in the Born approximation and found the expressions 
of the energy-loss rates for all the processes above. 



Later on, Brinkmann and Turner 19 calculated the axion emission rates 



in the nondegenerate (ND) limit and for a general degeneracy, for all three 
processes (a -c). They also could check the result of Iwamoto for the equal- 
nucleon cases in the (D) limit. For the NN interaction they used, however, 
constant nuclear matrix elements. The same results for the ND limit was 



also obtained previously by Turner in [18J. They reached the conclusion 
that ND regime is a better approximation of the axion emissivities for the 
conditions characteristic for a newly born neutron star. Other calculations 



considering also a OPE potential where performed by Raffelt and Seckel [20 
They studied the axion emission rates of iViV — > NNa processes in order 
to determine their Sa{&) structure function for NN interactions in neutron 
stars. With their calculations they concluded that the inclusion of pion mass 
effects do not reduce the axion emissivities by more than 50% even for ND 
regime. Our results will show that the contributions of pion mass to the 
coresponding emissivities due to nuclear effects are temperature dependent, 
and for a certain temperature interval are larger than this, as we will see 
later. Therefore, in this article we pay special attention to the ND regime 
but we present also the results for D regime. For the contribution of the np 
process we consider different chemical potentials for neutrons and protons, 
and as a result we can span different degeneracy degrees for the two species. 
Analogously to the case of neutrino pair emission (j2l|- 26j), one of the 



main difficulties for the calculation of axion emissivities is the appropriate 
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treatment of the strong NN interaction. In Iwamoto's calculations for the D 
limit the effects were included by replacing the nucleon momenta by their 
Fermi values in the angular part of the phase-space integrals. 

In the present work we include the full dependence on nuclear momenta of 
the nuclear matrix elements (NME) in the calculations of the axion emission 
rates by the NN bremsstrahlung processes (a - c). Our results for e a AT7v sep- 
arate explicitly the part corresponding to constant NME, which corresponds 
to the high-momentum limit of the previous works 12]- 22|). from the 



part including the missing nuclear effects due to the nucleon momenta depen- 
dence of the NME to facilitate the comparisons and further applications. We 
compare our results with those obtained by Brinkmann and Turner ([l9j|) 
for the ND limit, while for the degenerate regime we compare with those 
obtained in Ref. (Q) and ((l9|). 

2. Calculations 

The axion emission rate by NN bremsstrahlung is given by Fermi's Golden 



Rule formula (see for instance 191 ]) 



taNN = (271-)' 



n}- 



d 3 Pa 



(2vr) 3 2^J (2vr) 3 2£ a 



E a (Sx S|M| 2 ) 5\P)F(f) (1) 



where F(f) = fifg(l — fs)(l — fi) is the product of Fermi-Dirac distribution 
functions of the initial (1,2) and final (3,4) nucleons, fi = ^exp +1 

In Eq. (1) pj and E^ (i=l,4) are the nucleon momenta and energies, while 
p a and E a are the corresponding axion quantities; S is a symmetry factor 
taking into account the identity of the particles (1/4 for nn and pp channels 
and 1 for the np channel) and /ij are the chemical potentials of the nucleons. 

2 

In the non-relativistic limit E{ ~ m + j*-. Using the non-dimensional 
quantities 



19 



y = fi/T (fi = ji — m) and Ui = pj/2mT, the expressions 
of the Fermi-Dirac functions read fi = (exp Ui ~ Vi ■ The degenerate (D) 

limit satisfies y » 1, while in the non-degenerate (ND) limit y << —1. For 
S x S|M| 2 we use the following expressions for the nuclear matrix elements 
(in the OPEP approximation): 



Sxm\ 2 = Sx™-gW 



f 



— ■ M NN 



(2) 
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where 

9ai = Cm N /(f a /N) 
C is a dimensionless factor of order unity, which is model dependent, 



m n ~ m p = 940MeV = m N -> g an ~ g ap = g ai = C ■ 5.64 • 10 



-10 



(f a = 10 10 GeV is the Peccei-Quinn symmetry breaking scale, N = 6 repre- 
sents the color anomaly of the Peccei-Quinn symmetry). 
For the nn and pp the momentum-dependent factors Mnn read 

M ( |k|2 V | ( " |2 V I C-^W-I'I' (3) 
^|k| 2 + m;/ ^|l| 2 + m;y (|k| 2 + mj)(|l| 2 + mj) V ' 

with (5 = 3((k • l) 2 ) (k, 1 being the coresponding unit vectors for k and 1), 
while for the np process 

( Ikl 2 V / 111 2 \ 2 lkl 2 -lll 2 

y|k| +m^y y |1| +m^y (|k| +m*)(|l| + m;) 

where k = p x — p 3 and 1 = px — p 4 are the nucleon direct and exchange 
transfer momenta, respectively. The last (exchange) terms in the above ex- 
pressions arise from interference of two different reaction amplitudes. They 
contain contributions from the scalar product (k • l) 2 , which have been esti- 



mated ([12J,[19|]) by replacing them by their average values (denoted by (3) 
in the phase-space. There are two numerical values for f3 in the literature: 
(3 = 1.0845, in Ref. and (3 = 1.3078 in Ref. 0, but this difference 



(explained by Raffelt and Seckel) produces changes in our final results of 



only 1%. Since we compare our results especially with those of Ref. [19 
we use that value for (3. Thus, from kinematical constraints f3 =0 in the D 
regime, while it is 1.0845 in the ND regime (see the expression above). We 
have used for the np process, the NME of jToJ, with equal coupling constants 
for protons and neutrons. 



We follow the procedure of Brinkmann and Turner [19j to derive the ND 



limit, by performing the transformation to the center-of-mass system 

Pi + P2 Pi - P2 

p+ = — 7\ — ; p = — n — ; Psc = pa - p+; P4 C = P4 - p+ 
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=^Pi = p+ + p -; P2 = p+-p -; Ps = p+ + P3 C ; P4 = p+ + P4 C (5) 

From these relations and the conservation of momentum (axion momentum 
is neglected) we find p 4c = — p 3c . We define also the dimensionless quantities 



Ui = -^L(i = 1,4); u + = -5±-; w_ = — ; u 3c = (6) 

P+P P+P3c P P3c 

COS71 = 1 r; r; cos7 c = - r; r; COS7 = 1 r; (7) 

|P+I|P-| |P+l|P3c| |P-||P3 C | 

From the definition of the u variables above, and the conservation of 
energy, one can easily deduce the following relations 

«i,2 = u+ + u- ± 2(u + m_) 1/2 cos7i; u 3A = u + + m 3c ± 2(u + m 3c ) 1/2 cos7 c ; 

u_ = U 3c + ^a/2T (8) 

Let us now address the OPE potential. Following the method used in our 
previous papers ( j23| , [24j] ) , and after some lengthy algebra, one can express 
the matrix element M nn (eq.(4)) in terms of the scalar combinations |k| 2 + 11| 2 

1 2 2 

and k • 1 . Finally we expressed these NME in the following compact form 



S x S|M| 2 = [(3 - /3) - \M nn \lud} 



(9) 



where 



a — r . r 2 

\M I 2 ~m 2 * (10) 

\ Ivl nn\ nu d ~ ^ Q _ JJ . (J2 _ £j . (J4 \ > 

The coefficients A nn , B nn , C, D and E of eq.(10) are polynomials de- 
pending on the parameters m, T and m n and of variables U- and m 3c (for 
their full expressions, see Appendix A). Also we used the notation C$ = 
cos^icos^c + siwyxsiwYcCostfi, with (ft the angle between the vectors p + and 
P- 

Thus, the contribution of NME is split into a constant term,obtained also 
by Brinkmann and Turner ([l9|), Raffelt and Seckel(j2o|])) - which represents 
just its high-momentum limit (i.e. the limit to which the expression (3) con- 
verges when the pion mass is neglected compared to the nucleon momentum 
transfer) and a reduction term (see Appendix A) to be evaluated. After an 
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approximation which is numerically accurate within 1%, we succeeded to in- 
tegrate the expression of |M nn | 2 over the angles and finally we could express 
the axion emission rate in the ND limit in the following form 



rND l (NN) 

t aNN — t aNNy v ) I 1 ~ /g _ 

where 



ND _ ND (n\ I -■ *nucl\-- ) I fll\ 
<t\ A a \ S • I \ (3 — f3)I ND ' 



e™ N (0) = 2.68 x lO-y at eVT^(f/m n y (12) 



is the expression calculated by Brinkmann and Turner [19[ and Iq and 
InJd(NN) are double integrals over w_ and u 3c 



OO /■« 



Io D = / / y/(v-U3c)(u- - u 3c ) 2 e- 2u -du.du 3c (13) 
Jo Jo 



irm 2 r °° ru ~ 



JO 



V (2w_ + 2u 3c + m x ) 2 1 

with N = n, p, mi = m 2 /mT. 

A similar procedure for the np process yields 



S x S|M| 2 = 2 J^!k (i_y [ (7 _ 2/3) _ lMnp \l ud ] 



(15) 



where 



and ^ a 7v = [(7 — 2/3)/3]#^ is the effective axion nucleon coupling for the np 
case (see expression A.l in ref . (l9| ) . 

and final expressions analogous to Eq.(ll) and (12): 

ND JVD/ n \ ( 1 Inucl( n P) 



ND e M(Q) I -nuclV^J ^js 

c anp c anpl w ; ^ x (7 — 2(3)1$° J 

6^(0) = 2.68 x W-y aN ey^m 2 - 5 T^(f/m n r (18) 
The correction integral of the third term I^ t (NN) is replaced by 
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nm 2 poo r u- 



(2w_ + 2u 3c + mx) 2 J K J 

Using the same procedure, we calculated the emissivities for the nn(pp) 
processes in the D limit. In this case, in performing the integrals over 
the momenta and energies, we used the method of integration adopted in 



Refs.[9|,[17( and [25(. The expressions of the emissivities are 

e° aNN = & N {0) (l - (20) 



where 



or 



£ ?™(0) = ^rn^m-/gUV (21) 

(in natural units and Pf{N) = 2mTy) 
with N = n or p, and 

jD = f S ,:, f ' ^ j (22) 



n , s Kml f 00 f x f o/4tt 2 + 2 2 \ {6(2x + z)+7m 1 \ , , 

I°JNN) = £ / / z 3 ( — — x -) — ; J dxdz 

nuclK } mT J J Q \ e *-l) \(2x + z + m 1 ) 2 J 

(23) 

Similar calculations for the np process yields a final expression analogous 
toEq. (17): 



77~ D / 



with 

31^2 



C(°) = ^^ 2 - 5 ^ 6 - 5 ^ 4 ^/ 4 ^ /2 (l - Ay/2y m ) (25) 



3780tt 
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(expression obtained in ref. [191]). where y m = (yi+y 2 )/2 } Ay = \yx — 7/2 1/2, 
and the correction integral 



2 ,00 / 47r 2 + z 



mT 



</o 



x 



/14(2, + ,) + 17mA ^ (26) 



With these expressions at hand, we discuss the results in next section. 



3. Results 

In Figs. 1 and 2 we plotted the dependence on temperature of the relative 
corrections of the emissivities, Ae/eo = (eo — e)/e (where eo is the emissivity 
calculated in the high momentum limit - see Ref. [19J, while e is the emissiv- 
ity determined with our method), for ND regime (Fig. 1 - nn/pp processes 
and Fig. 2 - np process) due to pion mass effects (see Eqs. (11), (17)). One 
observes that for this regime, the relative corrections to the previous results 
for the emissivities are quite important and temperature dependent. For all 
processes, the corrected emissivities are reduced with 30% to 85%, depend- 
ing on temperature. Figs. 3 and 4 show the same relative emissivities, this 
time for the D regime (Fig. 3 - nn/pp processes and Fig. 4 - np process). 
In this case, the corrected emissivities are reduced at most with 11%. So, 
for this regime the contributions from pion mass effects (see Eqs. (20), (24)) 
do not seriously affect the corresponding emission rates (a result qualita- 
tively reported in [2(|). The dependence of these corrections on temperature 
is smoother in this case than for the ND regime. For more confidence, we 
give in Table 1 the absolute values for Ae/eo for ND and D regimes. We 
remark a very similar behaviour of the nn/pp and np relative corrections, for 
both regimes. 



9 



J [ivi e v \ 


/A,/, \nn/pp 


( \ c l c \ n P 
(/\e/e ) ND 


l [ivi e v \ 


( A c 1, \nn/pp 


(Ae/e ) D 


or 
ZO 


U.OOUo 


n q a on 

u.o4yu 


1 


U.lUoU 


n 1 c\'yc\ 
U.1U (U 


30 


0.7189 


0.7180 


2 


0.0960 


0.0997 


35 


0.6234 


0.6227 


3 


0.0900 


0.0931 


40 


0.5503 


0.5497 


4 


0.0846 


0.0874 


45 


0.4923 


0.4918 


5 


0.0800 


0.0824 


50 


0.4455 


0.4450 


6 


0.0756 


0.0780 


55 


0.4070 


0.4066 


7 


0.0720 


0.0740 


60 


0.3738 


0.3735 


8 


0.0686 


0.0705 


65 


0.3459 


0.3455 


9 


0.0653 


0.0673 


70 


0.3224 


0.3220 


10 


0.0627 


0.0643 


75 


0.3010 


0.3008 


20 


0.0440 


0.0435 



Table 1: Relative correction to the emissivities (relative emissivities) for ND and D regime, 
for all NN bremsstrahlung processes (nn,pp,np). Here Ae/eo = (eo — e )/ £ o, with e - the 
emission rates determined with our method and eo - the emissivities previously obtained 
inRef. 



(Ae/s) N 




1 1 1 1 1 1 1 1 1 1 1 1 

20 30 40 50 60 70 80 

T[MeV] 



Figure 1: Relative correction to the emissivities (Relative emissivities) due to pion mass 
effects for the ND regime, nn,pp processes. 
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Figure 2: Relative emissivities due to pion mass effects for the ND regime, np process. 
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0.07 
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(Ae/e)™ 



nn/pp 



— I — 
10 



— I — 
15 



— I — 

20 



T[MeV] 



Figure 3: Relative correction to the emissivities due to pion mass effects for the D regime, 
nn,pp processes. 
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Figure 4: Relative emissivities due to pion mass effects for the ND regime, np process. 



4. Conclusions 

In this paper we presented a method of determination of the NN axion- 
bremsstrahlung emissivities, for all processes, in both ND and D regimes, 
based on the inclusion of the full dependence of the NME on nuclear mo- 
menta. 

Starting with the ND case, worked out for the nn, pp and np cases, 
we found substantial reductions to the emissivities due to the combined 
effects of momentum-dependence in the temperature-dependent integrals - 
Ens. (14), (19). While these effects were somewhat foreseen in previous works 
([l9j), here we present an explicit calculation and quantitative results. We 
claim, based on the present results, that all axion mass limits that employed 
e a NN should be revised (the ones used in numerical simulation supernova 
codes), especially those for which the ND limit is highly relevant, such as 
SN1987A neutrinos ((l2l|). As it stands, the suppression of e a 7VAr can be im- 
portant and this feature postpones one to obtain a firm bound on m a for 
both popular KSVZ (0,0) and DSVZ (0,@) axion models. 

For the D case, our method allowed to check that the leading terms (Eq. 
(21), (25)) coincide with Iwamoto (jj|,[l7j]) and Brinkmann and Turner fflil]) 



results (see for example Eqs.(2.6) and (2.8) of Ref. [17]]), and also that the 
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reductions are up to 11% for all three processes (a-c). 

We also mention that other nuclear effects like spin-density fluctuations 
([27]], [2^]) or short range effects (TPEP) can further reduce the coresponding 
axion emissivities.We are currently calculating the NN axion-bremsstrahlung 
emissivities by including two pion exchange effects (through a TPEP) that 
can be mimicked by the exchange of a p meson ([29]). These effects are im- 
portant at distances below 2 fm ([3(|). Preliminary results show us a further 
reduction of the axion emission rates compared with the case presented here. 
The calculations will be reported in a future paper. For the moment it is fair 
to point out that apparently minor sources of error become actually impor- 
tant for the problem, already at the OPE approximation level, as explicitly 
shown above. Also, it is worth to mention that our method might be used 
to improve calculations for other physical processes of neutrino and axion 
emission in NS, in which strong interactions are also involved. 
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A. Appendix A 

We present here the correction matrix elements for all the cases and 
processes, which have not been presented in section 2 (Calculations). Starting 
with ND regime, for nn and pp processes, the correction matrix element is 



A nn = 2(3 - f5)m\Ul + 5(3 - (3)mlm 2 t U 2 + + 4(3 - (3)m 4 n m t U + (A2) 



M, 



r |2 
nn I nucl 




(Al) 



where 



B nn = 4m t 2 C/ p [2(3 - (3)m t U + + (7 - /3)m*] 
C = m\U\ + km\m\\J% + §m\m 2 t Xj\ + Amlm t U + + m\ 



'7T 



8 



(A3) 
(A4) 
(A5) 
(A6) 



D = 8m 2 t U p (m t U + + ml) 2 
E = 16miU p 
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with 

m t = 2mT; U + = u_ + u 3c ; U p = u_u 3c (A7) 

In ND case, for the np process we obtained the correction matrix elements 
and the corresponding coefficients (C, D and E are the same) as follows: 

■ , 2 _ 2 A np — C np ■ BnpCff, — D np ■ 

\Mnp\nucl ~ _ D . C 2_ E . C 4 + Q _ J) . £2 _ E . (J4 K A -*) 

where 

A np = (7 - 2(3)[2m 3 t Ul + 5m> t 2 ^ + m^m 2 + Am t U + )] (A9) 

B np = \2m 2 t U + Ul /2 {m t U + + ml) {A. 10) 

C np = 4[2(7 - 2p)mlU + U p + (17 - 2/3)m t 2 m2[/ p ] ( A n) 

D np = 48m t 3 f/ p 3 / 2 (A 12) 

For the D regime, the expressions for the correction matrix elements and 
for the corresponding coefficients are obtained by taking (3 = in previous 
relations (A.1-A.12). 
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